Ocular dominance patterns and the wire length minimization: a 
numerical study 
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We study a mathematical model for ocular dominance patterns (ODPs) in primary visual cortex. 
This model is based on the premise that ODP is an adaptation to minimize the length of intra- 
cortical wiring. Thus we attempt to understand the existing ODPs by solving a wire length 
minimization problem. We divide all the neurons into two classes: left- and right-eye dominated. 
We find that segregation of neurons into monocular regions reduces wire length if the number of 
connections to the neurons of the same class (intraocular) differs from the number of interocular 
connections. The shape of the regions depends on the relative fraction of neurons in the two 
classes. We find that if both classes are almost equally represented, the optimal ODP consists 
of interdigitating stripes. If one class is less numerous than the other, the optimal ODP consists 
of patches of the less abundant class surrounded by the neurons of the other class. We predict 
that the transition from stripes to patches occurs when the fraction of neurons dominated by the 
underrepresented eye is about 40%. This prediction agrees with the data in macaque and Cebus 
monkeys. We also study the dependence of the periodicity of ODP on the parameters of our 
model. 
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I. INTRODUCTION 

In the primary visual area (VI) of many mammals, 
most neurons respond to the stimulation of two eyes un- 
evenly: they are either right or left eye dominated. In 
some species, right/left eye dominated neurons are segre- 
gated and form a system of alternating monocular regions 
known as the ocular dominance pattern (ODP) (Wiesel 
and Hubel, 1965, 1969). In others, ODP is not observed 
(see Horton and Hocking, 1996b for a comprehensive list 
of species). ODPs, when observed, vary significantly be- 
tween different species and even between different parts 
of VI in the same animal. 

Most modeling studies of ODP (Erwin et al., 1995; 
Swindale, 1996) have addressed its development. They 
succeeded in generating ODPs of realistic appearance. 
However, several why rather than how questions remained 
unanswered. For instance, (1) why, from functional point 
of view, do the ODPs exist? (2) Why do some mam- 
malian species exhibit ODPs while others do not (Horton 
and Hocking, 1996b; Livingstone, 1996)? (3) Why do the 
monocular regions have different appearance (stripes as 
opposed to patches) between different parts of VI within 
the same animal (LeVay et al., 1985)? 

The question of functional significance of ODPs has 
been addressed theoretically using the wiring economy 
principle (Mitchison, 1991; Chklovskii, 2000). The idea 
is that evolutionary pressure to keep the brain volume 
to a minimum requires making the wiring (axons and 
dendrites) as short as possible, while maintaining neu- 
ronal functional properties (Cajal, 1995; Allman and 
Kaas, 1974; Cowey, 1979; Cherniak, 1992; Young, 1992; 
Chklovskii et al., 2001; Koulakov and Chklovskii, 2001). 
In many cases these functional properties are specified by 



the rules of establishing connections between neurons, 
or wiring rules. The problem presented by the wiring 
economy principle is therefore to find, for given wiring 
rules, the spatial neuronal layout that minimizes the total 
connection length. This approach allows to understand 
many features in cortical maps, such as orientation pref- 
erence maps (Koulakov and Chklovskii, 2001), as evo- 
lutionary adaptations, which minimize the total cortical 
volume. 

The goal of this study is to find the simplest model, 
which on one hand is supported by experimental evi- 
dence, and on the other encompasses most of OD phe- 
nomenology. The use of the simple model allows us to 
explore its parameter space completely and to give an- 
swers to the set of questions above. We also evaluate the 
dependence of the ODP period on the parameters of our 
model and compare it to the ODP periodicity observed 
in macaque monkey. We find that the experimentally ob- 
served variation of the period is in agreement with the 
wiring economy theory. 



II. MODEL AND METHODS 

A. Description of the model 

For the purposes of minimizing the cortical wiring we 
consider only intra-cortical connections since they con- 
stitute the majority of gray matter wiring (LeVay and 
Gilbert, 1976; Peters and Payne, 1993; Ahmed et al., 
1994). We therefore disregard the thalamic afferents and 
other extra-cortical projections. In an attempt to make 
wiring economy argument more quantitative, we propose 
a model describing the component of intracortical cir- 
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cuitry sensitive to OD. The principal elements of our 
model are therefore the connection rules between corti- 
cal neurons. To assess the sensitivity of the intracortical 
wiring to OD we examine the connections in the cortical 
layer AC/3, where OD is most strongly pronounced. Such 
sensitivity has been studied by Katz et al. (1989). They 
made three observations regarding the wiring rules: 

i) Neurons in the layer AC (3 near the interface between 
two OD columns arborize more in home rather than in 
the opposite eye column. Therefore neurons establish 
more connections with the neurons dominated by the 
same rather than by the opposite eye. 

ii) Axons and dendrites of these neurons have a ten- 
dency to bend away from the interface between OD 
columns. This implies that not only they avoid penetra- 
tion to the opposite OD column but also they attempt to 
maintain sufficiently high number of connections in the 
home column. 

in) Axons or dendrites penetrating through the oppo- 
site eye column to the next same eye column are never 
observed. This means that retinotopy has little effect 
on connections in layer AC/3. Indeed the neurons on the 
edges of two same eye columns separated by one opposite 
eye column have on average receptive fields centered next 
to each other. If connections in AC/3 were sensitive to the 
retinotopic coordinates, these two edges should be con- 
nected (Mitchison, 1991). However out of 21 cells exam- 
ined Katz et al. (1989) observed none producing axons 
reaching the next same eye domain. The only possibility 
for such cells to be connected is due to the overlap be- 
tween dendritic and axonic arbors of two cells separated 
by more than 500 /im. Such possibility is small because 
of the strong repulsion of the connections by the opposite 
eye column located between two cells (observation if). 

These three observations lay the basis of our model 
which we now describe. The elementary unit of our model 
mimics the columnar organization of the cortex (Mount- 
castle, 1957) and uniformity of ODP along the direction 
normal to the slab. The elementary unit is therefore a 
microcolumn, which is defined as a box, spanning the 
cortex perpendicular to its surface, whose other two di- 
mensions are smaller than the characteristic scale of ODP 
(ss 500/i), and yet large enough to include many neu- 
rons. A possible choice of dimensions for such a micro- 
column is thickness of cortex (~ 1.5mm ) x 30/i x30/i, 
in which case it includes about 310 cells in VI (Rockel et 
al., 1980). The microcolumn units are therefore arranged 
on a square lattice with 30/i period. 

Although the choice of the elementary unit size may 
seem arbitrary, the results of our calculation are inde- 
pendent of the choice. The size of the unit is analogous 
to the integration step, which does not affect the value 
of an integral significantly if chosen to be small enough. 

Motivated by the second observation in layer AC/3 
listed above, i.e. that neurons maintain a fixed number of 
connections in the home OD column, we make the follow- 
ing assumption about the connection rules. Each micro- 
column unit must establish connections with N s distinct 
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FIG. 1 Our model illustrated. The units dominated by the 
left and right eyes are shown by the full and empty circles 
respectively. Each unit is required to make N 3 connections to 
the units of the same OD and N connections to the units 
of the opposite OD. In this illustration N s = 4 and N — 2. 
The connections satisfying these rules are shown for two units, 
right and left eye dominated. Small numbers of connections 
are chosen for the ease of illustration; in actual implemen- 
tation of the model both N s and N are large (see below). 



units dominated by the same eye and N a units domi- 
nated by the opposite eye. These connection rules are 
illustrated in Fig. ^ 

Only the relative values of N s and N Q (rather than 
absolute) are important because of the arbitrariness in 
the definition of microcolumn. Thus, if a 30 x 30/i mi- 
crocolumn receives N s — 10 4 projections from the same 
OD column and _/V = 10 3 projections from the oppo- 
site OD column, a 60 x 60/i microcolumn receives four 
times less projections respectively, or N s — 2.5 • 10 3 and 
N Q = 2.5 • 10 2 . This is because with coarser units each 
projection is more effective: connecting to one 60 x 60/i 
unit implies connecting to four 30 x 30/i units. Both 
these implementations of the model produce the same 
OD pattern, discretized in "pixels" of different size. The 
important quantity, which is invariant with respect to the 
change of "pixel" /microcolumn size, is the ratio between 
N a and N Q (equal to 10 in this example). This is the first 
parameter of our model. 

The second parameter of our model is the filling frac- 
tion of the units (microcolumns) dominated by the left 
eye afferents Jl with respect to the total number of 
units, averaged over several ODP periods. This param- 
eter is the fraction of the left eye dominated units 
by fjj + fji = 1. For the majority of important cases 
fh — frt. — 1/2, however on the periphery of visual field 
one of the eyes (ipsilateral) is underrepresented. There- 
fore, its filling fraction is less than 1/2. 

The third observation above implies that the compo- 
nent of connections sensitive to OD is not sensitive to the 
retinotopy, and both numbers N s and N do not depend 
on the position of the receptive field of the unit. This may 
be due to significant scatter of the receptive field of the 
cells within cortical column on the scales of about 1mm 
(Hubel and Wiesel, 1974). The position of the receptive 
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field of the microcolumn is therefore vaguely defined and 
cannot affect OD sensitive connections significantly. 



B. Methods 

Given these wiring rules we look for an optimal lay- 
out of the microcolumn units which minimizes the total 
length of connections. To find the layout minimizing the 
total wirelength we employ a combination of computa- 
tional and analytical techniques. To make our choice of 
methods clear we first comment on the expected proper- 
ties of the solution. 

A possible solution of our model is the Salt and Pep- 
per layout in which the units dominated by right and left 
eyes are uniformly intermixed. In this layout the units 
belonging to different eyes are not segregated, ODP is 
not formed, and the local values of the filling fraction are 
equal to 1/2 (by local value we understand an average 
over a domain including many units yet small compared 
to the period of ODP). It should be contrasted to the 
case when units dominated by the same eye fill in large 
domains i.e. form the ODP. In the latter case the local 
values of the filling factor of each eye vary from to 1. 
However one can imagine an intermediate situation when 
the local filling fraction varies from 1/2 — a to 1/2 + a, 
where the amplitude of variation < a <C 1/2. This 
corresponds to the case of weak segregation into ODP. 
The weak segregation is found in squirrel monkey where 
ODP has fuzzy appearance and until recently was sus- 
pected not to be formed (Horton and Hocking, 1996). If 
a = 1/2, i.e. the local filling fraction varies from to 
1, the ODP's have sharp appearance. Using the general 
terminology from binary mixtures (cite diblock copolimer 
paper) we call this regime the strong segregation limit. 

The methods useful in the strong segregation limit are 
not good in the weak segregation case and vice versa. We 
use the simulated annealing to find the optimum phases 
for the strong and nearly strong segregation cases. Hav- 
ing found the optimum phase in the strong segregation 
case to assess the period of ODP we use the exact enu- 
meration technique, which compares layouts belonging to 
the same phase with different periods. The treatment of 
the weak segregation case requires the use of continuous 
variables and is done employing the perturbation theory. 
Below we describe these methods in more detail. 



1. Simulated annealing 

The parameters of Metropolis Monte-Carlo method 
(Metropolis et al., 1953) are optimized to render most 
consistent results for multiple restarts. We use square 
20 x 20 array of units with periodic boundary conditions. 
The units are either left or right eye dominated. At each 
step the algorithm attempts to change the dominance of 
one unit to the opposite. The value of the average filling 
fraction /^n is enforced by adding the following term to 



the total connection length: 
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where L and are the current values of the total wire- 
length and average filling fraction. Such term in the func- 
tional keeps the current value /j, close to the required 
value f La . 

To map out the phase diagram the values of Jl change 
from 0.2 to 0.8 in 0.02 increments. The values of N s and 
N D satisfy the condition N s + N Q = 30 and are changed 
in unit increments, i.e. have the following values: 12, 18; 
13, 17; 14, 16; 15, 15; 16, 14; etc. The phases at the 
intermediate points are taken from the nearest points, 
where result is available. 

The Monte-Carlo temperature is gradually annealed 
from 0.24L/A to 0.008L/A (N = 400 is the total number 
of units) in 5000 sweeps through the entire system (20 x 
20 x 5000 steps). The resulting layout is then examined 
and the phases visually identified. 



2. Perturbation theory 

Salt and Pepper layout is relatively easy to study due 
to its uniformity, and can be solved exactly (Chklovskii 
and Koulakov, 2000). If a layout does not deviate signif- 
icantly from Salt and Pepper, i.e. the weak segregation 
case takes place, it can also be solved exactly. This im- 
plies that the wire length can be written as an explicit 
functional of density distribution of the units. Such func- 
tional was evaluated and optimized with respect to the 
density variations by Chklovskii and Koulakov (2000). 
The optimization shows that ODPs are formed for the 
values of parameter | A^ s — N a \/N s > 0.02. However the 
simulated annealing method cannot distinguish weak seg- 
regated ODP from Salt and Pepper for \N S — N Q \/N S < 
0.2. There are two reasons for the failure of simulated 
annealing to do so: 

• Simulated annealing is performed at small but fi- 
nite temperature that destroys weakly segregated 
ODP. 

• The units can be either completely right or left eye 
dominated. This implies that OD can change only 
sharply in the described annealing version. This is 
useful for obtaining the strongly segregated phases, 
which occupy major part of the parameter space. 
However, in the weak segregation limit the local OD 
changes smoothly. Thus used version of simulated 
annealing performs poorly at N s w N a . 

We therefore replace the simulated annealing results 
by those from Chklovskii and Koulakov (2000) at small 
values of parameter N s — N Q (see the phase diagram be- 
low). 



4 



3. Calculation of the ODP period. 

To evaluate the period of ODP precisely, we first deter- 
mine the phase (Salt and Pepper, Stripes, or Patches) for 
the given set of parameters N S /N Q and fx,, using meth- 
ods described above. We then take a lattice containing 
a large number of units, which exceeds sufficiently the 
lattice used in simulated annealing. This is possible be- 
cause the method of determining period described below 
is much less time consuming than simulated annealing. 
We then arrange the two types of units on the lattice, 
using ODP determined by the simulated annealing, and 
vary the period of the pattern to find the period produc- 
ing the minimum of the wire length. Below we describe 
the procedure for both Stripes and Patches in more de- 
tails. 

i) Stripes 

To find the optimum period for stripes we use an array 
containing 300 by 300 units. This array includes three pe- 
riods of the stripes, which run parallel to one of the sides 
of the region. Each period therefore includes 100 units, 
containing nx left and tir right eye units, nx+nxt = 100. 
By varying nx we accomplish the change in the filling 
fraction of the left (ipsilateral) eye, according to the for- 
mula: fx=nx/(nx+nxt.)- We consider a string of 100 
units at the center of the array, which is representative of 
all the units in the configuration. For each of the central 
units the computer program establishes connections, ac- 
cording to the connection rules. Most of the calculations 
are done for N S + N Q = 300. We check that results change 
for different N s + N Q in a predictable fashion (see be- 
low, Results, subsection IIII.D.TJl . Stripes therefore have 
a fixed period in terms of number of units (100). To 
find the optimal spatial period of the stripes we vary the 
shape of each elementary cell in the 300 by 300 array. 
Thus, if the rectangular cell dimensions are a x perpen- 
dicular and a y parallel to the stripes, we vary both a x 
and a y , keeping the area of elementary cell a x a y = 1 
constant. By doing so we do not change the density of 
units, but vary the spatial OD period, according to the 
formula A = lOOa^. For each value of a x the cells are 
reconnected according to the connection rules. Special 
care is taken about exclusion of the boundary effects by 
making sure that none of the units on the edges of the 
array is connected to. After the optimum period is found 
the period in terms of number of units is changed from 
100 to another value, closer to the value of spatial period, 
to check for the absence of geometric artifacts, associated 
with distortions of elementary cells. The change of the 
spatial period after this procedure is typically absent but 
in extreme cases does not exceed 3%. 

ii) Patches 

Since our results indicate that a triangular crystal of 
Patches is formed (see Fig. EP), we consider an array in 
the shape of parallelogram commensurate with the tri- 
angular arrangement of Patches. The lattice sites in the 
array, representing units, are also arranged on a triangu- 
lar lattice. Their positions are given by x(i,j) = i + j/2 



and y(i,j) — jV3/2, where i and j are integers varying 
between 1 and 51. Here I is the period of ODP to be op- 
timized. The centers of Patches are located at points 
x c n,m = In + lm/2 and y c (n,m) — lm^/i/2. Each 
patch includes lattice sites at the distance from a cen- 
ter determined by the filling fraction of the ipsilateral 



eye: R = ly fxV3/2Tt. The units within/outside the 
patch are left/right eye dominated. The units in the 
configuration are then represented by the central par- 
allelogram: i,j = (21 + 1)..3/. For each of the units con- 
nections are made according to the connection rules with 
N s + N a = 300. The optimum period is obtained by 
varying parameter I. 



4. Fourier analysis of the ocular dominance patterns 

To determine the experimental dependence of the ODP 
period on the filling factor, the image of ODP in macaque 
monkey (Horton and Hocking, 1996a) is converted into 
a digital format. In this format the image is represented 
by a set of pixels. A pixel with coordinates x and y is 
represented by a number s(x,y), equal to for the right 
eye dominated and 1 for the left eye dominated area. For 
each position in the image we then determine the local 
value of the average filling fraction of the ipsilateral eye 
and the value of local OD period. Both these calculations 
are similar. To do the calculation at a certain point in 
the map, given by coordinates (xo,yo), we surround the 
corresponding pixel by a square, containing 64 x 64 pixels 
(black square in Fig. [21 3.7 x 3.7mm). The dimensions of 
the square are such that one hand it contains a few ODP 
periods (about 3), which is needed for averaging, and on 
the other hand it is small compared to the characteristic 
dimensions over which the properties of ODP change (~ 
lem, see Fig. |5J). To determine the filling fraction we 
average the scanned image over the square: for position 
(xo,yo) in the map the local value of the average filling 
fraction is given by 
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To determine the local value of ODP period we perform 
the Fourier transform of the s (x, y) — fx in the square. As 
a result we obtain a set of numbers s(q x ,q y ), representing 
the Fourier transform amplitudes, defined on a 64 x 64 set 
of wave vectors (q x ,q y ). The spectral power, represented 
by \s (q x , q y )\ 2 , is shown in Fig. [3] for one of the points 
in the pattern, corresponding to Stripes. It clearly has a 
bimodal appearance, indicating the average in the square 
direction of the stripes. We then determine the average 
value of the wave vector, using the formula: 



(<l(xo,yo)} 



(3) 
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a) 



FIG. 2 The image of the striate cortex of macaque monkey 
1 left hemisphere from Horton and Hocking (1996a). The 
left /right eye dominated areas are shown by black/white. For 
the pixels in the shaded area we evaluate the filling factor and 
OD period, displayed in Fig. 1171 below. The dashed square 
gives an example of the region containing 64 by 64 pixels, for 
which the filling fraction and Fourier transform are calculated. 
It has dimensions 3.7 by 3.7 mm. 
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FIG. 3 The spectral power for a point in the pattern occu- 
pied by stripes. The spectrum has the bimodal appearance, 
characteristic of stripes. The direction of the modes is deter- 
mined by the direction in which ODP changes (perpendicular 
to stripes). The distance of the modes from the center deter- 
mines the local value of ODP period by Eq.Q. The spectral 
power in the scale bar is in arbitrary units. 



The value of the mean ODP period is then defined as 
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b) 




This value for each pixel in the shaded area in Fig. [21 is 
shown in Fig. IT7I 



FIG. 4 Ocular dominance patterns for ft = 1/2 and N a = 
N — 4. (a) A realization of the Salt and Pepper phase gives 
minimal wire length (/ ~ 9.67 lattice constants per neuron), 
(b) A realization of the Stripe phase is suboptimal (I w 10.24). 



III. RESULTS 

A. Small number of connections 

We start by finding optimal layouts for three illustra- 
tive examples of wiring rules with small numbers of con- 
nections, N s and N a . We caution the reader that because 
of the small numbers of connections phase assignments 
may seem arbitrary. These examples are chosen to illus- 
trate our main results which will be confirmed both in the 
lattice model with large N s and N a later in this section 
and in the continuous model (Chklovskii and Koulakov, 
2001). 

For the first two examples we set equal numbers of left 
and right dominated neurons, /& = fa = 1/2. In the first 
example each neuron connects with equal numbers of the 
same-eye and other-eye neurons, N s = N Q — 4. Then the 
optimal layout is the "chess board" of left/right neurons, 
FigQp,. This layout is a realization of the Salt and Pepper 
phase, FigU^L, because each neuron has an equal number 
of left and right neurons among its immediate neighbors. 
To calculate the length of connections per neuron, I, we 
notice that in this layout all neurons have the same pat- 
tern of connections. By considering one of them, Fig^, 
we find that I = 4 + 4^2 w 9.67. This layout is optimal 
because each neuron makes all of its connections with 
immediate neighbors. 

A suboptimal layout for the same wiring rules is illus- 
trated by a realization of the Stripe phase, Fig^D. In 
this layout each neuron has the same pattern of connec- 
tions up to a mirror reflection. By considering one of 
them, Fig0J}, we find 1 = 6 + 3V2 » 10.24, greater than 
I ~ 9.67 for the Salt and Pepper phase. Here each neuron 
has among its immediate neighbors only three other-eye 
neurons, while the wiring rules require connecting with 
four other-eye neurons. A connection to a more distant 
neighbor is longer making the layout suboptimal. We 
confirm the optimality of the Salt and Pepper phase for 
N s = N for large N s , N both numerically and analyti- 
cally. 

In the second example each neuron connects with more 
same-eye than other-eye neurons: N s = 5, N = 3. Then 
a realization of the Salt and Pepper phase, Fig[5Ji is not 
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a) 




b) 




FIG. 5 Ocular dominance patterns for ft, — 1/2 and N s — 5, 
N = 3. (a) A realization of the Salt and Pepper is subopti- 
mal (I fa 10.24). (b) A realization of the Stripe phase gives 
minimal wire length (I w 9.67). 
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FIG. 6 Ocular dominance patterns for ft, — 1/4 and N s — 5, 
N = 3. Realizations of the (a) Salt and Pepper (I w 11.26) 
and (b) Stripes (I w 11.49) are suboptimal. (c) A realization 
of the L-Patch phase gives minimal wire length (I ~ 10.67). 



optimal anymore. The length of connections per neuron 
is I w 10.24, while the Stripe phase, Fig|SjD gives I « 
9.67. The Salt and Pepper phase loses in wiring efficiency 
because there are not enough same-eye neurons among 
immediate neighbors and connections with more distant 
neighbors are needed. The Stripe phase, FigEb rectifies 
this inefficiency by having each neuron make connections 
only with immediate neighbors. Thus, clustering of same- 
eye neurons is advantageous if each neuron connects more 
with the same-eye than with the other-eye neurons. 

In the third example we use the same wiring rules 
(N s = 5, N Q = 3) but take different numbers of left/right 
neurons, /l = 1/4, fa — 3/4. The realizations of the 
Salt and Pepper phase is shown in Fig[fjL and of the 
Stripe phase in FigHJb- In these layouts, different neurons 
have different patterns of connections. To find the wiring 
length per neuron we average over different patterns and 
find for the Salt and Pepper phase I ~ 11.26 and for the 
Stripe phase I « 11.49. A more efficient layout is the L- 



Patch phase, Fig[U, where I w 10.67. Although we can- 
not prove that the L-Patch phase is optimal, this seems 
likely. Thus, the optimal shape of monocular regions de- 
pends on the relative numbers of left/right neurons. 



B. The shape of OD columns 

After giving some examples of ODPs with small num- 
bers of connections N s and N a we discuss the opposite 
case of large numbers. As we show below in Section 
IIII.D.ll the shape of OD columns in this case does not 
depend on the absolute values of parameters N s and N a . 
It is determined by the ratio N s /N and by the rela- 
tive amount of ipsilateral neurons fi (assuming that the 
left eye is ipsilateral). Depending on the values of pa- 
rameters N S /N Q , and /l, optimal layout belongs to the 
one of the eight phases shown in Fig. [7J where ipsilat- 
eral and contralateral-eye dominated neurons are shown 
by black and white regions respectively. These phases 
can be divided into three major classes. The first class is 
represented by the unsegregated Salt and Pepper layout, 
in which two types of neurons are uniformly intermixed 
(Figure [7J\) . The second class includes Stripy layouts, 
shown in Figures [7p, E, G, I. The third class consists of 
Patchy layouts, displayed in Figures 03, F, H, G. 

We distinguish several subclasses of Stripy phases. 
First, it is the sharp Stripes (Figure[7J), which consists of 
alternating lamellar monocular regions. Second, it is the 
weakly segregated Stripes (Figure EP,E). In this ODP 
the variation of density of left /right eye dominated neu- 
rons is small. This is an intermediate pattern between 
the unsegregated Salt an Pepper and the sharp Stripe 
layouts. This phase is therefore fragile and difficult to 
obtain numerically. In some cases simulated annealing 
can produce such a phase, Figure [TJH. In the other cases 
the weak segregated phase can only be obtained by the 
perturbation theory, which can carefully account for a 
weak variation of neuronal density. Such case is shown 
in Figure [7p. Third, we also obtain Stripy phases that 
show a tendency to become Patches, by e.g. their longi- 
tudinal modulation, such as shown in Figure Ep. 

Similar subclasses exist among Patchy layouts. We 
obtain sharp, weakly segregated (obtained from simu- 
lated annealing or perturbation theory), and elongated 
Patches, which are shown in Figures 0J, F, D, and H 
respectively. Finally, there are mixed phases containing 
both Stripes and Patches, such as in Figure [7)3. These 
ODP's arc shown on the phase diagram (PD) in Figure|Sl 
The phase diagram shows the optimum phase (minimiz- 
ing the total wire length) for given values of parameters 
N s /N and f L . 

The important feature of the PD is its left-right eye 
symmetry. It is apparent from the symmetry of Figure |H| 
with respect to the line fi = 1/2. This is a consequence 
of the left-right eye symmetry of our model, implying 
that the connection rules, defined by numbers N s and 
N a are independent on whether a neuron is left or right- 
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FIG. 7 The phases obtained by perturbation theory (C and 
D) and by simulated annealing (the rest). Each simulation 
array of 20 by 20 units in reproduced 4 times in each Fig- 
ure. The color bar is the key for Figure |S] A: Salt and Pep- 
per, B: Stripes mixed with Patches; C and D: weakly segre- 
gated Stripes and Patches obtained by the perturbation the- 
ory; E and F: weakly segregated Stripes and Patches obtained 
by simulated annealing; G: modulated Stripes; H: elongated 
Patches; I: sharp Stripes; J: sharp Patches; 



eye dominated. For this reason the phase for > 1/2 
can be obtained from the point with the same N s /N a 
and the value of the filling fraction equal to 1 — /x < 
1/2. This corresponds to the replacement of the white 
regions in Figure by black and vice versa. A similar 
correspondence takes place when one compares ODP's 
in left and right hemisphere. 

Another important feature of the PD is the existence of 
the Salt and Pepper region around the line N S /N Q = 1. 
This implies that the difference between N s and N a is 
the driving force of segregation into ODP. The larger the 
difference, the sharper the ODP becomes. 

The area of the PD adjacent to fj, = fn = 1/2 is 
occupied by Stripy ODPs. At small values of the filling 
fraction the phases are Patchy. A transition from Stripes 
to Patches occurs at /z, « 0.38 almost independently 




16/14 
15/15 
14/16 
13/17 
12/18 



FIG. 8 The phase diagram showing the optimum phase for 
given values of parameters /l and N a /N . For the color key 
see Figure [7] 
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FIG. 9 The example of transition from Patches to Stripes 
at N s /N — 3. A, B, and C show the optimum phases for 
Jl = 0.34, 0.36, and 0.38 respectively. 



on parameter N s /N a . An example of such transition for 
N s /N a = 3 is shown in Figure 03 

The reasons for the transition for small values of 
N S /N D — 1 are discussed in Koulakov and Chklovskii, 
1999. For /arge?differences, when ODP becomes sharp, 
the transition occurs due to the presence of surface con- 
tribution to the wire length. To demonstrate this we 
present the following argument, which is rigorously valid 
in the asymptotic limit of large number of connections to 
the same-eye neurons, i.e. N s » N . In this limit con- 
nections to the same-eye neurons are the most abundant 
and therefore most costly, from wire length prospective. 
Hence, we can disregard connections to the opposite- 
eye neurons in the first approximation. Consider a unit 
near the interface between two OD columns (top unit in 
Fig. II Op t . The connection arbor of this unit to the same 
OD units, shown by empty circles in Fig. ^3 is strongly 
biased toward the home column, since the unit has to 
maintain certain number of connections there. This effect 
has been seen by Katz et al., 1989, in macaque striate cor- 
tex (see also the discussion in the Model Section above) . 
For units away from the interface the connection arbor 
within the same OD column is close to a circle (Fig. ITU1 
bottom unit). Thus, circular arbor renders the minimum 
wirelength in the absence of constraints, such as the in- 
terface between OD columns. With the interface present 
the connection arbor to the same eye neurons is therefore 
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FIG. 10 Connection arbors of two units in the array: near 
the interface of OD columns (top) and in the center of a col- 
umn (bottom). The units dominated by left and right eyes 
are shown by squares and dots respectively. Two units for 
which the arbors are displayed are shown by stars. They are 
connected to other units, which are encircled. The parameters 
for this particular layout are: N 3 = 23, N = 7, ft = 0.40. 
The connection arbor of the boundary unit (top) to the units 
of the same OD is significantly deformed, compared to the 
corresponding connection arbor of the unit in the middle of 
the column (bottom). This is similar to the observation of 
Katz et al., (1989). This deformation gives rise to the sur- 
face cost associated with formation of the interface between 
columns. 



strongly deformed with respect to the optimum. Hence, 
the presence of the interface between the OD columns 
implies an increase in the wirelength, and is therefore as- 
sociated with a surface cost (similar to Malsburg, 1979). 
This surface cost drives the transition from Stripes to 
Patches. Indeed if /l *C 1/2 the patchy phases have 
much shorter length of the surface compared to Stripes. 
This is because Patches shrink when /j, — ► reducing 
their surface length, whereas Stripes become narrower, 
keeping their surface length the same. However, this 
is not true for fi, — fn = 1/2 where Stripes have a 
shorter surface for numerical reasons. Therefore, when 
jh is decreased, the Stripes should unavoidably condense 
into Patches to minimize the surface cost. This is similar 
to droplets of water assuming circular shape to minimize 
the surface energy. 

We conclude therefore that two important transitions 
occur in our model. 

• The transition from unsegregated Salt and Pepper 
to weakly segregated and then sharp ODP is driven 
by the difference between parameters N s and N a 
characterizing the intra-cortical circuitry. 



• The transition from Stripy to Patchy ODP is driven 
by the decreasing filling fraction of the ipsilateral 
eye and occurs at fj, w 0.4, if left eye is underrep- 
rcscntcd. 



C. Transition from Stripes to Patches: comparison to 
experiments 

Our phase diagram in Fig.|S|shows that the transition 
from Stripes to Patches occurs when fi, ss 0.4 for a wide 
range of N s /N . This value will be compared now with 
the value of /l at which the transition occurs in the ex- 
periment, found from the relative area occupied by left 
eye dominated neurons. The conclusion that the Patch 
phase becomes optimal when contralateral eye dominates 
is, indeed, non-trivial, because there may be a system of 
alternating wide and narrow monocular stripes instead. 

We test our conclusion on the data from macaque mon- 
key first (Horton and Hocking, 1996a). The relative area 
occupied by the left /right eye depends on the location 
in VI. In the parafoveal part of VI both eyes are rep- 
resented equally, i.e. Jl ~ 0.5. ODP has a stripy ap- 
pearance, in agreement with the phase diagram. Away 
from the foveal region contralateral eye becomes domi- 
nant. The ODP becomes patchy there (LeVay et. al., 
1985), just as expected from the theoretical phase dia- 
gram. We verify the location of the transition by using 
the following algorithm. We find fi for each point of the 
pattern by calculating the relative area occupied by the 
left /right regions in a window centered on that point and 
a few OD periods wide (dashed lines in Fig. Ill|) . Then 
we draw a contour corresponding to — 0.4, Fig. ^2 
We observe in Fig. ^] that stripes indeed become patchy 
at the black contour indicating /r = 0.4. 

In Cebus monkey the ODP has a similar transition 
(Rosa et al., 1992). For monkey C06L from Rosa et 
al., 1992, we determine visually that along the horizon- 
tal meridian the transition occurs at the eccentricity of 
20 — 40 deg. According to the plot of the relative repre- 
sentations given in Rosa et al., 1992, /l changes in the 
range 0.32 — 0.42 at these eccentricities. Our theoretical 
conclusion about a transition at /l = 0.4 falls into this 
interval. For the upper 45 degree meridian of the same 
monkey the transition occurs at the eccentricity of 30—40 
degrees or at filling fractions 0.33 — 0.43. Again, the pre- 
dicted value belongs to this interval. We conclude that 
these data are consistent with the results of our model. 

In cats the ODPs have a patchy appearance (Anderson 
et al., 1988; Jones et al., 1991). In this case our theory 
implies that one of the eyes should dominate. According 
to some reports (Shatz and Stryker, 1978; Crier et al., 
1998) the filling fraction of the contralateral eye in cat 
VI is about 0.8 in young animals (before P22). This may 
lead to Patches in cat VI. The strong contralateral bias 
disappears in older animals (Crier et al., 1998). This is 
consistent with other reports (Anderson et al., 1988) that 
both eyes are represented almost equally in adult cats. 
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FIG. 11 Transition between the Stripes and Patches occurs 
at theoretically predicted value /t. Shown is a fragment of 
the macaque ODP from Horton and Hocking (1996a). Areas 
dominated by the left eye are grey and neurons dominated 
by the left eye are white. The black contour corresponds 
to Jl = 0.4 averaged over a window, whose dimensions are 
shown by the dashed square (3.7 x 3.7mm). The points of 
transition from Stripes to Patches coincides with the black 
contour. 



D. The period of ocular dominance pattern 

1. Scalability of the Model 

One of the general features of our model is that the pe- 
riod of OD pattern becomes larger, when the total num- 
ber of connections is increased. Indeed, the size of the 
connection arbors grows if both N s and N a are increased, 
given that the density of units (l/(30/i x 30/i)) is kept 
constant. This is because one has to go further to find 
the necessary number of connections to satisfy the wiring 
rules. Since dimensions of the connection arbors set up a 
characteristic scale for the OD pattern, the period of the 
latter also increases. This property of our model, which 
we call scalability, is discussed in this subsection. 

Let us define scalability in a mathematically rigorous 
fashion. Assume that one has found a minimum wire 
length configuration for certain set of parameters N s , N Q , 
and /x (or fu = 1 — Jl). Assume that both N s and N a 
are very large. In this case the dimensions of connec- 
tion arbors are much larger than the lattice spacing, and 
one can ignore the fine structure imposed on the con- 
nection arbors by the square lattice. This is exactly the 
limit in which our model has some validity, both because 
realistic numbers of neuronal connections are large and 
because we would like to avoid artifacts introduced by 
the square lattice. An example of connection arbors for 



a neuron for some set of N s and N a are shown in Fig. IT21 
(left). This neuron and its connection arbors resemble 
the top neuron, marked by the star, in Fig. The 
connection arbors in Fig. ll2l look like continuous circular 
pieces, due to the large N s and N Q limit (square lattice 
makes the boundaries of connection arbors look like stair- 
cases, whose steps are too small to show in the picture 
for large N s and N a ). Imagine now a geometric transfor- 
mation, in which the dimensions of the connection arbors 
of all of the neurons, as well as the OD pattern itself, are 
blown up by the same scaling factor rj > 1. The two- 
dimensional density of the neurons must be preserved 
during this transformation. The obtained new OD pat- 
tern and the new connection arbors are shown schemat- 
ically in Fig. El (right). Since the transformation is ap- 
plied to the two-dimensional objects, and each of the di- 
mensions was stretched by the factor rj, each neurons in 
the new configuration will receive rj 2 N s and rj 2 N a con- 
nections from the same and opposite eye neurons. The 
newly obtained configuration (Fig. ^| right) will satisfy 
wiring rules with connection numbers given by t] 2 N s and 
rj 2 N . Note that the filling fraction is not changed by 
this transformation. It is /l = /_r = 1/2 in Fig.^J Will 
this be the minimum wire length configuration for this 
set of parameters? 

To prove that the new configuration minimizes the to- 
tal wire length for the new set of parameters rj 2 N s and 
i] 2 N we notice that the total wirelength for the new 
configuration is given by rfL, where L is the total wire- 
length for the original configuration (Fig.ll2llefth This is 
because each neuron now receives the number of connec- 
tions increased by rj 2 , and each connection is stretched 
by rj. Imagine now that one finds a new configuration, 
which has the total connection length L' < rfL. Let us 
take this more optimal configuration and shrink it by the 
factor of rj. We obtain a configuration, satisfying wiring 
rules for the set N s and N a , whose total wirelength is 
L'/rj 3 < L. But this contradicts to our postulate that 
the original configuration in Fig. 1121 (left) is optimal for 
the set of parameters N s and N a . Thus the stretched 
configuration provides the minimum of the wirelength 
for the new set of parameters rj 2 N s and rj 2 N . 

This property is important, because once the solution 
for given N s and N a is found, one can obtain other so- 
lutions, with the set of parameters r] 2 N s and 7] 2 N Q , by 
stretching the original configuration by the factor of r/ 
uniformly in all the directions. The important property 
which remains the same for all these related configura- 
tions is the ratio between the numbers of the same and 
other eye connections, N S /N Q . Thus, we conclude that 
this ratio determines the shapes of the OD patterns, 
which is unchanged during the uniform stretching pro- 
cedure. 

What is changed in the uniform stretching is the ODP 
period? Since the period is proportional to the stretch- 
ing parameter ry, while the total number of connections 
is proportional to rj 2 , we conclude that the period is pro- 
portional to the square root of the total number of con- 
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FIG. 12 The original ODP and connection arbor for one 
of the neurons, marked by the white dot (left). Compare to 
Fig. 1101 top neuron. The stretched configuration is shown on 
the right. 



nections, when the ratio N s /N is kept constant, 
implies that 

A(N s ,N ,f L )=D-X(N a /N ,f L ) 



This 



(5) 



Here D = 2a^J{N s + N a ) /n ~ r], where a = 30/i is the 
size of the microcolumn unit. Parameter D has a mean- 
ing of the average axonal arbor diameter. The quantity 
X(N S /N , /i) is the normalized OD period, calculated in 
the units of the average axonal diameter. This quantity 
is introduced here for easier comparison to the experi- 
ment. Notice that this quantity does not depend on the 
total number of connections. The latter dependence is 
entirely absorbed by the parameter D. 

Scalability is valid for the limit of large N s and N a , 
when square lattice effects can be ignored, and our model 
becomes continuous. Does scalability apply to our nu- 
merical solution, described in subsection lII.B.3t ' To check 
this we plot the ratio A (N s , N a , f£) /D, obtained using 
methods described in lII.B.3l for different values of the to- 
tal number of connections N s + N in Fig. El If Eq. JHJ is 
valid, this ratio should not depend on the total number of 
connections. As evident from Fig. 1131 this property is in- 
deed satisfied. Hence, below in this subsection we always 
present the results for X(N S /N , fa) = A/D, which are 
independent on the total number of connections, assum- 
ing that the latter dependence can be easily recovered 
using Eq. (J5J. 



2. Dependence on parameter N s /N 

We now examine the dependence of normalized period 
X(N S /N a , fi) [see Eq. (JSJ] on the parameter N S /N Q , for 
Jl = 1/2, when we have to consider the stripe phase, 
according to subsection IIII.BI The results are shown in 
Fig. [21 These results have been obtained using methods 
described in subsection lII.B.3l In general, the normalized 
period increases with increasing parameter N S /N Q . This 
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FIG. 13 The independence of the ratio A/D [see Eq.@] on 
the total number of connections. This implies that the lattice 
effects, violating scalability, are insignificant. The data are 
obtained for /l = 1/2. Circles, triangles, and diamonds show 
results for N s /N equal to 4/3, 17/3, and 14 respectively. 



increase in the OD period can be understood consider- 
ing the interplay between connections to the same and 
opposite eye units. Indeed, the presence of connections 
between the same eye units implies affinity between the 
same OD neurons. An increase in the relative number 
of such connections (N s /N ) strengthens such affinity. 
The OD columns provide a neighborhood rich with the 
same eye neurons; even more so, on average, for coarser 
domains. Thus stronger affinity between the same eye 
neurons (N S /N D ) leads to an increase in the size of OD 
domains. This effect is produced by wiring economy prin- 
ciple, since the latter provides a substrate for the affinity 
of connected neurons. 

We now examine Fig. ^] in more detail. The relative 
period diverges in the limit N s >> N Q . The divergence 
can be described by the asymptotic formula 



X(N s /N ,f L 



(6) 



shown in Fig. 1141 by the dotted curve. The divergence 
can be understood as follows. Imagine that the neurons 
do not have to connect to the neurons of the opposite 
OD, i.e. parameter N Q — 0, N s ^ 0, which corresponds 
to the extreme case N s >> N Q . In this case the opti- 
mum wire length configuration consists of only two large 
domains, dominated by left and right eye neurons, occu- 
pying a half of VI each. This is because of the notion 
of surface contribution, introduced in subsection IIII.BI 
To minimize this interface contribution the system phase 
segregates into two large domains. Thus, in the case 
N Q = ODP has maximum possible period, spanning 
the entire VI. This explains the tendency of the period 
diverge in the limit N a ^ {N s /N = 00) in Fig. [21 
What happens if N Q 7^ 0? Since the neurons now have 
to connect to the opposite eye neurons, phase segregated 
configuration (two large domains spanning the entire VI) 
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is no longer optimum. The system introduces more inter- 
faces between OD columns to shorten distances between 
opposite eye neurons. More interfaces implies a reduc- 
tion in the OD period. Thus, finite N s /N a leads to the 
finite OD period. This is reflected by the asymptotic 
dependence © and the dotted curve in Fig. [21 

An interesting phenomenon observed in Fig. 1141 is the 
presence of a few discontinuous changes in the OD pe- 
riod. One such a change occurs at N s /N a w 1.15 and 
consists in an abrupt increase in the OD period by a 
factor of about 2.3. Another discontinuous transition oc- 
curs at N s /N a i=s 1.20 and the corresponding increase in 
the period is by a factor of 2. Note that these transi- 
tions are truly discontinuous, i.e. they occur at discrete 
points. They appear smooth in Fig ll4l due to a sparse 
sampling (the real data points are shown by dots) . Note 
also that the quantity D in Eq. © changes negligibly 
in the interval between N S /N Q = 1.1 and 1.2, which im- 
plies that both OD period A and the normalized period 
A change approximately by the same factor. Such dis- 
continuous changes in the OD period in the weakly seg- 
regated regime, i.e. when the ODP is not well defined, 
may be responsible for the coarsening of ODP in strabis- 
mic squirrel monkeys (see Discussion for more details). 



3. Dependence on the filling fraction 

The dependence of the normalized period 
X(N S /N , Jl) on the filling fraction of the left eye 
Jl is shown in Fig.^| These results have been obtained 
using methods described in subsection III. B. 31 Four 
dependencies are shown, for four values of the parameter 
N 8 /N marked on each curve. The general tendency for 
the period to grow with increasing parameter N a /N 0l 
described in the previous subsection, is evident in the 
Figure. 

For small values of parameter N s /N the period in- 
creases when the filling fraction moves away from fa = 
1/2, as long as one stays within the same phase (Stripes 
or Patches). This behavior is seen for the two bottom 
curves in Fig. 1151 In the transitional region between 
Stripes and Patches the OD period experiences a discon- 
tinuity, marked by the dotted lines. For large values of 
N S /N D , the dependence of the period on Jl is opposite: 
the period decreases, as the filling factor deviates from 
1/2, as demonstrated by the top curve in Fig. 1151 

We now compare this behavior of our model to the 
observations in the striate cortex of macaque monkey 
(Horton and Hocking, 1996a), using Fourier transform 
method described in subsection III. B. 41 To make this 
comparison possible the following technical consideration 
is taken into account. The Fourier transform method 
evaluates the average value of the spatial frequency of 
the ODP (Q). The OD period is then calculated by the 
formula A = 2ir/ (Q). For the Stripe phase this proce- 
dure results in the value, which is close to the period of 
stripes. For Patches it results in the distance between 




FIG. 14 Dependence of the normalized OD period, de- 
fined in Eq. on the parameter N a /N . The data points 
are shown by dots connected by lines. The dotted curve 
shows the asymptotic fit obtained for the values of param- 
eters N s /N >> 1 [Eq.©]. Two horizontal dotted lines show 
the plateau values of the period separated by discontinuous 
transitions at N s /N ~ 1.15 and 1.20. 
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FIG. 15 The normalized OD period A as a function of the 
filling fraction of one of the eyes. Four dependencies are shown 
for the values of parameter N a /N equal to 2, 3, 5, and 11 
correspondingly. The central segment of the dependencies, 
between /l = 0.4 and 0.6 corresponds to the Stripe phase. 
The dependencies in the regions between 0.2 < /l < 0.35 
and 0.65 < fn < 0.8 have been calculated for Patches, as 
indicated in the Figure. In the regions of transition between 
Stripes and Pathes the curves are connected by dotted lines. 



rows of patches, which is smaller than the period by the 
factor V3/2 » 0.87 (see Fig. |TBJ|. Thus, to make com- 
parison to the experiment possible, the values in Fig. 1151 
corresponding to Patches should be multiplied by the fac- 
tor 0.87. The result of this procedure is shown in Fig. 1171 
by the gray line. 

Fig. El shows that the period observed in the exper- 
iment decreases when the filling factor of the ipsilateral 
eye deviates from 1/2. This warrants the use of the top 
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FIG. 16 For Patches Fourier transform produces distance 
between rows, rather then the period of the OD pattern. The 
distance between rows is a natural successor of the period of 
Stripes after the transition to Patches occurs (see also Fig. [§J . 
The distinction between Patch period and distance between 
rows should be taken into account for accurate comparison to 
the experimental observations. 
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FIG. 17 Comparison of the OD period observed in the exper- 
iment (Horton and Hocking, 1996a) in macaque striate cortex 
(dots) to the results of our theory (gray curve). The former 
is obtained using Fourier transform method described in sub- 
section The latter is the top curve in Fig. 1151 with the 
sector of the data corresponding to Patches corrected by the 
factor v3/2 ~ 0.87 for compatibility with the Fourier trans- 
form. The only fitting parameter used is D — 0.46mm [see 
Eq. ©]. 



curve in Fig.Elto represent the theoretical result. Since, 
the shape of the theoretical dependence does not change 
much when N S /N Q > 10, the parameter N s /N a cannot 
be established from the comparison of the theory to the 
experiment. To obtain the gray curve in Fig. 1171 the 
N 8 /N = 11 dependence in Fig. 1151 was multiplied by the 
fitting parameter D = 0.46mm. This is the only fitting 
parameter used. As seen in Fig. El °ur theory describes 
the experimentally observed dependence quite well. 

The widths of the ipsilateral and contralateral eye 
stripes in macaques has also been measured indepen- 
dently by Tychsen and Burkhalter (1997). Based on their 
results one can evaluate the ODP period and the filling 
fraction: 




ipsi 



FIG. 18 The dependence of the ODP period on the filling 
fraction based on data from Tychsen and Burkhalter (1997). 
The results for non-strabismic adult macaque monkeys are 
presented by markers: monkeys M25 right VI (diamonds), 
M25 left VI (circles), 418 (squares), and 906 (dots). The 
curve shows the best parabolic fit (see text). 



Here fi ps i, Wj, and Wc are the filling fraction of the 
ipsilateral eye, and the ipsilateral/contralateral column 
widths respectively. The dependence of the period on 
the filling fraction can therefore be established. This de- 
pendence is shown in Fisr llSI 

The best parabolic fit to the data in FigEDis given by: 



A(/)=A(l/2) l-a(/, 



1/2)' 



(8) 



A = Wi + Wc , h 



Wi/ {Wi + W c ) . (7) 



The coefficient a = — 6.94±6.38 is estimated using boot- 
strap (Efron and Tibshirani, 1993). The expectation 
value of the coefficient is therefore below zero, as seen 
from FigEO The probability of the coefficient to be 
greater than zero is p — 0.11 as evaluated by bootstrap 
procedure, which is used since the distribution of coeffi- 
cients a is non-gaussian. This implies that with great de- 
gree of certainty one can assume that the period of ODP 
does decrease with the filling fraction deviating from 1/2. 

It should be noted that the value of coefficient a can be 
obscured by the variability of ODP period between differ- 
ent individuals, since data in FigEDare assembled from 
three monkeys (four VI 's). To reduce the impact of inter- 
individual variability we then normalized the period for 
each area VI by the mean value for each individual ani- 
mal. The value of the coefficient is then a = — 7.98±6.24, 
with the probability of positive coefficient p — 0.055. 
Thus the decrease of the period with filling fraction is 
even more certain, when the inter-individual variability 
is accounted for. The value of coefficient a obtained from 
the theory is 2 (Fig. El N s /N a = 11). It is consistent 
with both measurements. 
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IV. DISCUSSION 

In this work we propose a model which can account for 
most of experimentally observed features of ODPs. Our 
model has two principal parameters. The first parame- 
ter characterizes the intracortical circuitry. It is the dif- 
ference between the number of connections to the same 
and to the opposite OD neurons. Our results suggest 
that this difference is the driving force of segregation 
into ODPs. We argue therefore that the sensitivity of 
the intra-cortical connectivity to OD provides a reason 
to formation of OD columns (see below). The second 
parameter is the fraction of neurons dominated by the 
ipsilateral eye. This parameter determines the shape of 
monocular regions in ODP. In the majority of the pri- 
mary visual area of macaque and Cebus monkeys this 
parameter is close to 50%, which implies that both ipsi- 
and contralateral eyes are equally represented. However, 
in the proximity of monocular crescent the ipsilateral eye 
becomes underrepresented. This is because the inputs 
into the eye are blocked by the nose of the animal. Our 
theory suggests that the decrease in the filling fraction 
of the ipsilateral eye drives the transition in the ODP 
structure from stripy (zebra skin like) to patchy (similar 
to leopard skin). The transition occurs when the fraction 
of the ipsilateral eye dominated neurons approaches 40% 
in both macaque and Cebus monkeys (see below). We 
also analyze the dependence of OD period on the param- 
eters of our model and find satisfactory agreement with 
experimental data. 



A. On the functional significance of OD columns 

Each neuron in our model establishes certain number 
of intra-cortical connections with neurons dominated by 
the same and the opposite eye. As suggested by experi- 
mental studies in macaque striate cortex, neurons in layer 
4C/3 typically make more connections with neurons of the 
same OD (Katz et al., 1989). Thus, from wiring econ- 
omy prospective, connections with neurons of the same 
OD are more important than the opposite eye connec- 
tions. Therefore, it is advantageous to form OD columns, 
since they provide environment rich with the same OD 
neurons, which results in shortening connections to the 
same eye neurons. The wiring economy principle thus 
provides a natural reason for the existence of OD pat- 
terns, i.e. answers the first question in the program listed 
in the Introduction. 

Our model suggests that in primates with weakly de- 
fined OD columns, such as squirrel monkey (Hubel et 
al., 1976; Livingstone, 1996; Horton and Hocking, 1996b) 
and owl monkey (Livingstone, 1996), the difference be- 
tween these two components of intracortical connectivity 
is not large. Such difference may be increased in these 
animals by experimentally induced strabismus. Indeed, 
strabismus reduces correlated activity between opposite 
OD cortical neurons, therefore reducing their connectiv- 



ity N Q . Reduction in N Q unbalanced by the correspond- 
ing reduction in N s increases the parameter N s /N a and 
leads to sharpening of OD columns, according to our 
phase diagram in Fig. [S] Such sharpening is indeed ob- 
served experimentally (Shatz et al., 1977; Lowel, 1994; 
Livingstone, 1996). This phenomenon was also predicted 
theoretically by Goodhill (1993). 

The two parameters of intra-cortical circuitry, N s and 
N a , represent in our model the interplay between two 
classes of processing performed by the visual cortex. The 
first class includes the processing of the monocular im- 
age, for which connections to the same OD neurons are 
necessary. The second class includes various tasks related 
to stereopsis, which require comparison of two monocu- 
lar images, relying on the connections between the oppo- 
site OD neurons. We proposed above that the function 
of OD columns is to shorten the connections between 
the same eye neurons. The impact of stereoscopic vi- 
sion should therefore be the opposite: strong stereoscopy 
should make ODP weaker. Inversely, weak stereoscopy 
induces sharp ODPs. The latter statement is justified 
by the observations in animals with experimental stra- 
bismus. However, one should be careful about this state- 
ment, since the relation between functional and anatom- 
ical properties may not be direct. 

The situation in the animals with lateral eye posi- 
tioning, such as mice, rats, tree shrews, etc., is differ- 
ent. In such animals the visual pathway is almost com- 
pletely crossed, i.e. VI in each hemisphere is strongly 
dominated by the contralateral eye [Drager, 1974, 1975, 
1978; Drager and Olsen, 1980; Gordon and Stryker, 1996 
(mouse); Hubel, 1977 (rat); Casagrande and Harting, 
1975; Mully and Fitzpatrick, 1992 (tree shrew); Horton 
and Hocking, 1996b (other species)]. As suggested by 
Antonini et al. (1999) this implies that the ODP con- 
tains only two large monocular 'columns', each spanning 
the whole hemisphere. This can be interpreted as an 
OD having very large period, spanning both striate cor- 
tices. This picture can be fitted into the framework of 
our model. Indeed, we predict that if the number of con- 
nections to the other OD neurons (N a ) is very small the 
OD has a very large period fFiglTil). Thus our model 
predicts that the number of connections received by each 
neuron from the neurons of the same OD (N s ) is much 
larger than the number of opposite eye connections (N a ) 
in the species with lateral eye positioning. This should in- 
clude the cross-hemispheric projections. This statement 
should have functional consequences. Since N Q is small, 
synthesis of images from two eyes is weaker. Therefore 
the animals have to find another strategy to implement 
stereopsis. Hooded rats for example use vertical head 
movements to estimate distances (Legg and Lambert, 
1990). Our conclusion about small N a should also ap- 
ply to the superior colliculus in these animals (Colonnese 
and Constantine-Paton, 2001), as well as to the tectum in 
lower vertebrates (Schmidt and Ticman, 1985), in which 
cases the visual inputs cross over almost completely too. 

To summarize, our model encompasses most of the 
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phenomena related to the sharpness and observability of 
ODP. It relates the interspecies variability in the ODP 
to the relative amount of binocular interaction occurring 
in different species. Thus, we predict, that in the ani- 
mals with weakly segregated columns (squirrel monkey) 
N s ss N a ; in the animals with sharp columns (macaque) 
N s is much larger than N a [confirmed by Katz et al., 
(1989)]; and, finally, in the animals with lateral eye po- 
sitioning, N a should have a value, whose contribution to 
the wirelength can be neglected. 

B. Variation of the ODP period in the extrafoveal region 

Another consequence of a decreasing Jl in macaque is 
a decrease in the ODP period (LeVay et al., 1985). In 
Fig. El we compare the result of our theory to the data 
from macaque monkey (Fourier transform applied to data 
from Horton and Hocking, 1996a). From this comparison 
we conclude that, according to the wiring economy prin- 
ciple, parameter N S /N Q >> 1, or cells establish much 
more connection with the same OD cells, than with the 
opposite. This result of is consistent with the observa- 
tions of OD sensitive circuitry in the striate cortex of 
macaque by Katz et al. (1989). 

We chose the regions proximal to the horizontal merid- 
ian for this comparison. This is based on the assumption 
that cortical properties, such as N s and N , change lit- 
tle along this meridian. This assumption is in part sup- 
ported by the fact that OD periodicity changes little on 
the large segment of the meridian occupied by stripes, 
spanning the region between about 2 and 25 degree ec- 
centricity (notice very little scatter in Fig. 1 171 around the 
point /j, = 1/2). The changes in the period begin to 
occur when Jl deviates from 1/2. This is illustrated by 
Fig. El Other authors notice a decrease in the period 
when comparing vertical to horizontal meridian. Studies 
based of computer reconstructions report about 2 fold de- 
crease in OD periodicity comparing these areas (LeVay 
et. al., 1985), while others, based of flat-mounts (Horton 
and Hocking, 1996a), report a more moderate change. 
Such variation cannot be accounted for by a decrease in 
parameter fi, alone, since Jl is about 1/2 on both meridi- 
ans in close proximity to parafoveal region (j20 degrees of 
eccentricity). Our model suggests two possibilities based 
on the variation in the intracortical circuitry, described 
by N s and N Q . Since such differences in the circuitry 
may also be responsible for variability of the OD period 
between different animals, we discuss this possibility in 
the next subsection. 

C. Variability of the ODP period from individual to 
individual 

Studies in macaque monkeys (Horton and Hocking, 
1996a) reveal large inter-individual variability of the 
stripe period. The stripe period is 1072 ± 164/i along the 
VI border, after comparison of 6 animals. Two factors 



may contribute to this phenomenon in the framework of 
our model. (1) The basic diameter of axonal and den- 
dritic arbors D varies from animal to animal. This could 
be due to changes in N s , N a , or neuronal density. (2) 
The ratio between monocular and binocular interactions 
N s /N varies. The former reason is justified by Eq. JSJl. 
The latter can be understood from Fig. El Simply speak- 
ing, monocular interactions (N s ) favor formation of OD 
columns, making them wider, in an effort to provide same 
OD rich environment for all the neurons. Binocular inter- 
actions (-/V ) favor interfaces between columns, since in- 
terfaces bring opposite OD neurons closer to each other. 
They therefore decrease OD period. When N S /N Q in- 
creases the OD period increases too (Fig.EJ- This may 
occur when comparing different individuals. 



D. ODP period in strabismic animals 

Similar consideration may apply to the experiments 
in strabismic animals (Lowel, 1994; Livingstone, 1996). 
Since strabismus reduces correlations between eyes, its 
effect in our model is to reduce parameter N Q . Hence, 
the ratio N s /N is increased. According to our model 
(Fig. I14fl this generally leads to an increase in the rela- 
tive OD period (ratio of the basic OD periodicity to the 
connection range D). This result is understood from the 
interplay between affinity between the same eye neurons 
(iV s ), increasing the period, and the affinity between op- 
posite OD neurons (-/V Q ), reducing OD period. Since the 
latter is reduced by strabismus, the OD period grows. 

The degree of the period change depends on the de- 
crease in the number of interocular connection, and is dif- 
ficult to estimate. A curious feature displayed by OD pe- 
riod in our model is an abrupt increase at N S /N Q w 1.15 
by a factor of about 2.3, cf. Fig. El This implies that 
close to point N s /N a = 1.15 the OD period may be very 
sensitive to developmental manipulations. This finding 
may have correlate in squirrel monkey, for which the ob- 
served increase in period is by a factor 2.9 — 3.6 (Horton 
and Hocking, 1996b). These data are obtained from com- 
parison to a single strabismic animal. The following sce- 
nario is possible, comparing squirrel monkey to the stra- 
bismus experiments in owl monkey (Livingstone, 1996), 
in which no significant increase in periodicity is observed. 
Parameter N S /N Q in squirrel monkey passes the point 
N s /N a = 1.15 due to strabismus, leading to the substan- 
tial increase in period. In owl monkey parameter N S /N D 
is above 1.15 in normal animal. Strabismus therefore has 
little effect. This scenario is consistent with sharper OD 
columns in normal owl monkeys (N s /N a > 1.15) than in 
normal squirrel monkeys (N s /N a < 1.15) (Livingstone, 
1996; Horton and Hocking, 1996b). Experimentally in- 
duced strabismus in cat leads to an increase in the OD 
period by a factor of 1.3 (Lowel, 1994; Goodhill, 1993; 
see however Jones et al., 1996). Our model suggests that 
parameter N S /N Q > 1.15 in cat, and the increase in the 
period is due to the smooth part of the dependence in 
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Fig. 1141 which may not be so substantial as in squirrel 
monkey. 



E. On the importance of wiring minimization 

The relevance of wiring economy principle to the neu- 
ronal spatial organization can be illustrated by the fol- 
lowing thought experiment (Koulakov and Chklovskii, 
2001). Imagine taking a cortical area and scrambling 
neurons in that area, while keeping all the connections 
between neurons the same. Since the circuit is un- 
changed, the functional properties of the neurons remain 
intact. Therefore, from the functional point of view, the 
scrambled region is identical to the original one. In fact, 
the only difference caused by scrambling is in the length 
of neuronal connections. Therefore, it is hard (if not im- 
possible) to justify the existence of systematic cortical 
maps, such as OD pattern, without invoking the cost of 
making long neuronal connections. Although some the- 
ories of map formation may not explicitly mention the 
wiring optimization principle, it is present implicitly, usu- 
ally in requiring the locality of intra-cortical connections. 

How important is the constraint imposed by wiring 
minimization? In principle one can imagine development 
of an organism, which has 30% excess of wire with respect 
to the existing ones. It turns out that the existence of 
such an organism is close to impossible. Indeed, imagine 
that an external object, such as a blood vessel, is intro- 
duced in certain area of the gray matter. In this case 
some of the neuronal connection have to go around the 
vessel, therefore increasing in length. If the nerve pulses 
are to be delivered at the original speed and/or inten- 
sity, the elongated axons and dendrites have to be made 
thicker, to increase the pulse propagation speed and de- 
crease dendritic attenuation. This leads to more obsta- 
cles on the way of other neuronal connections and so on. 
Thus, introduction of a new blood vessel leads to an in- 
finite series of axonal and dendritic reconstructions. The 
same is true about the extra connection volume, resulting 
from wasteful neuronal positioning. Such infinite series of 
reconstructions can diverge, which implies that the con- 
nection volume resulting from more and more reconstruc- 
tions increases indefinitely. In this case the new blood 
vessel can never be inserted without sacrificing signifi- 
cantly the brain function. It turns out that mammalian 
brain has reached the verge of this so called wiring catas- 
trophe (Chklovskii and Stevens, 2001), so that it gets in- 
creasingly more difficult to accommodate excess volume 
in the nerve tissues. 

The wiring catastrophe occurs when the fraction of 
axons and dendrites in the cortical volume reaches 60%. 
Electron microscopy studies of cortical slices show that 
the actual volume occupied by neuronal processes is 
about 55% (Chklovskii and Stevens, 2001). The brain 
therefore has approached the barrier imposed by wiring 
catastrophe. Further increase of the volume fraction of 
neuronal processed may deteriorate the brain function. 



F. Comparison to other models 

As discussed in the previous subsection, wiring opti- 
mization is the only known way to relate neuronal lay- 
out (as manifested in the ODP) to the statistics of neu- 
ronal connectivity. Models of the ODP development that 
do not explicitly rely on wiring optimization invoke it 
implicitly, usually requiring the locality of intra-cortical 
connections. 

Here we discuss the relationship of our model to the 
models that invoke wiring optimization explicitly. 

In his pioneering work, Mitchison (1991) studied a 
question whether ODP minimize the wiring volume rela- 
tive to the Salt and Pepper layout. He assumed that the 
inter- neuronal connectivity is determined both by oc- 
ular dominance and retinotopy with all neurons having 
the same connectivity rules. He found that the answer to 
this question depends on the detailed assumptions about 
axonal branching rules. In particular it depends on the 
value of axonal branching exponent a. He has shown that 
if all axonal segments have the same caliber (a = oo), 
than ODP's are indeed advantageous for certain range of 
ratios of same-eye to opposite-eye connections. He also 
showed that if a < 4 than the ODP do not save wiring 
volume relative to Salt and Pepper. However, existing 
data seems to suggest that axonal caliber branches with 
a < 4 (Deschenes and Landry 1980, Adal and Barker, 
1965). 

The case of axonal branching with the cross-sectional 
area conservation corresponds closely to our model be- 
cause we require a separate connection for each neuron. 
The reason we find that ODP minimize wiring length 
is because we drop the retinotopy requirement on inter- 
neuronal connection rules, an assumption supported by 
the experimental data (Katz) . Although, effectively con- 
nections are roughly retinotopic, connection rules may 
vary from neuron to neuron thus providing some flexi- 
bility. The advantage of our approach is its simplicity 
allowing us to map out a complete phase diagram and 
make experimentally testable predictions. The full the- 
ory of the ODP will require a detailed analysis of axonal 
branching which must rely on better knowledge of axonal 
branching rules. 

Jones et al. (1991) proposed an explanation for why 
ODP have either Stripy or Patchy appearance. They 
assumed that neurons are already segregated into the 
ODP (by considering units whose size equals the width of 
monocular regions) and found that the difference between 
Stripy and Patchy appearances of the ODP could be due 
to the boundary conditions, i.e. different shape of VI in 
different species. Although the correlation between the 
shape of VI and the ODP layout is observed, the model 
of Jones et al. does not explain why peripheral repre- 
sentation of macaque VI has patchy ODP or why ocular 
dominance stripes run perpendicular to the long axis of 
VI in some parts of VI but not in others. Moreover, it is 
the local structure of ODP that is likely to determine the 
shape of VI and not the other way around. Therefore, 



16 



unlike Jones et al., our work proposes a unified theory 
of ODP including Salt and Pepper, Stripy and Patchy 
layouts, and is based on local inter-neuronal connectivity 
rules. 

Another model related to wiring length minimization is 
the elastic net model studied by Goodhill and coworkers 
(1993). The original formulation of the model minimized 
the cost function which penalized for placing nearby neu- 
rons whose activity is not correlated, a choice justified by 
computational convenience. Later the elastic model was 
generalized by the introduction of a C measure. Maxi- 
mization of C measure effectively corresponds to penal- 
izing for placing correlated neurons far apart. Unlike 
wiring optimization the penalty does not increase beyond 
a distance called cortical interaction. Because of this, 
elastic net often yields solutions where left and right eye 
neurons are completely segregated into left and right eye 
maps. 

Our wiring optimization models can be viewed as a 
sub-set of models described by C measure. The advan- 
tage of our wiring optimization approach is that it has a 
transparent biological justification for the cost of placing 
neurons far from each other - the cost of wiring. Because 
of this, wiring optimization is a natural choice for ques- 
tions related to the anatomy of intra-cortical connections. 

Wiring optimization provides a link between neuronal 
connectivity and spatial layout. Thus, it leaves open 
the connection between connectivity and computational 
function. Unlike most other models, which assume that 
neurons should be connected only if they are correlated, 
wiring optimization makes other assumptions about con- 
nectivity possible, for example connecting neurons with 
anti-correlated firing. 



G. Conclusion 

Our theory relates functional requirements on the neu- 
ronal circuits to its structural properties. In particular, 
our model relates the amounts of the neuronal intraoc- 
ular and interocular interactions, and the filling frac- 
tion of ipsilateral neurons, to the structure of the ocular 
dominance pattern. We conclude that the interspecies 
variability in the ocular dominance patterns may be ex- 
plained by differences in the underlying cortical circuitry. 
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